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SUMMARY 

The  objective  of  this  grant  was  to  investigate  the  problem  of  excitation 
of  a parallel-plate  EMP  simulator.  This  problem  has  been  studied  in  the  past 
by  a number  of  investigators  including  the  authors  of  this  report.  However, 
in  the  previous  investigations,  the  approach  taken  was  based  on  the  use  of 
the  leaky  mode  concept  wherein  the  fields  excited  in  an  open  waveguide  are 
represented  in  terms  of  the  leaky  modes  supported  by  the  guide.  It  is  well- 
known,  however,  that  the  .eaky  modes  are  not  proper  solutions  of  Maxwell's 
equations  for  the  entire  space,  and  that  in  order  to  be  complete,  the  leaky 
mode  representation  must  be  supplemented  by  the  contribution  from  the  contin- 
uous spectrum. 

The  motivation  of  the  present  effort  was  to  derive  a numerically 

il 

exact  solution  to  the  source  excitation  problem  with  a view  to  providing  a 
standard  for  comparing  the  approximate  leaky  wave  representation. 

Two  approaches  were  used  to  attack  the  source  excitation  problem  in 
the  open,  parallel-plate  waveguide.  The  first  of  these  is  described  in 
Attachments  A and  B and  is  based  upon  the  solution  of  a finite  Wiener- 
Hopf  equation  using  function-theoretic  techniques.  The  above  method  is 
most  useful  for  wide  plate  widths  and  is  in  fact  restricted  in  its  applica- 
tion to  this  case.  The  second  method  described  in  Attachment  C is  entirely 
numerical  in  nature  and  employs  the  method  of  moments  to  solve  the  integral 
equation  for  the  plate  current  induced  by  the  given  source.  This  approach 
complements  the  Wiener-Hopf  formulation  since  it  is  most  useful  for  small  to 
moderately  large  plate  widths. 

Sz 

All  of  the  results  given  here  are  for  an  assumed  variation  e in 
the  longitudinal  direction,  with  S,  the  wave  number  along  2,  as  the  important 
parameter.  Though  not  carried  out  here,  the  complete  solution  to  the 
field  excitation  problem  can  be  derived  by  performing  a spectral  integration 
with  respect  to  the  variable  3. 
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The  choice  of  Che  other  parameters,  e.g.,  plate  width  and  separation 
was  suggested  by  the  Air  Force  Weapons  Laboratory  and  the  LuTech  Company. 
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The  latter  organization  was  collaboratively  engaged  in  an  experimental  study 
on  simulators  with  the  research  group  at  Harvard.  During  the  course  of 
these  experiments,  certain  anomalous  behaviors  were  observed  for  which  no 
obvious  explanation  was  available.  It  is  hoped  that  the  results  presented 
here  would  shed  some  light  on  the  problem  of  understanding  some  of  the 
anomalous  bahavior. 

To  summarize,  this  report  presents  the  results  of  a study  of  the 
problem  of  source  excitation  of  an  open,  parallel-plate  waveguide  modeling 
the  EMP  simulator  and  presents  extensive  numerical  results  derived  by  using 
two  complementary  approaches  — one  for  the  small  width,  and  the  other  for 
the  wide-plate  case. 
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In  this  work  we  investigate  numerically  the  problem  of  the  source 
excitation  of  an  open,  parallel-plate  waveguide.  The  following  assump- 
tions are  made  for  the  source  current  1)  the  current  is  oriented  in  the 
y-direccion,  2)  it  is  located  at  x ■ 0,  3)  there  is  no  variation  in  the 
y-direccion,  4)  and  the  current  has  exp  (idz)  behavior  along  the  longi- 
tudinal z-direction.  We  provide  graphical  output  for  the  EM-field 
components  as  functions  of  a longitudinal  propagation  constant  and 
transverse  coordinates  and  then  discuss  these  results. 
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I.  INTRODUCTION 


In  the  previous  report  [1],  we  derived  analytical  expressions  for 
the  source  excitation  of  an  open  parallel-plate  waveguide.  However, 
these  formulas  were  very  complicated,  and  it  became  necessary  to  evaluate 
them  numerically.  The  purpose  of  this  report  is  to  present  the  numerical 
results.  The  computer  program  contained  in  Appendix  A was  written  and 
used  to  obtain  the  field  distribution  as  a function  of  the  longitudinal 
propagation  constant  and  the  transverse  coordinates.  The  numerical  out- 
puts are  presented  in  graphical  forms.  The  Cyber  175  at  the  University 
of  Illinois  was  used  for  all  of  the  numerical  studies. 

The  organization  of  the  report  is  as  follows:  Section  II  contains 
a statement  of  the  problem  and  the  basic  formulation.  Section  III 
presents  the  real  and  imaginarv  parts  and  the  amplitude  of  the  component 
field  distribution  as  functions  of  several  parameters  in  graphical  form 
and  a detailed  discussion  of  the  numerical  results.  Finally,  Section  IV 
is  the  conclusion. 
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II.  STATEMENT  OF  THE  PROBLEM  AND  BASIC  FORMULATION 


In  this  section  Che  fields  due  to  a vertical  current  located  inside 
an  open,  finite  waveguide  are  investigated.  The  geometry  of  the  problem 
considered  is  shown  in  Figure  1.  This  structure  consists  of  two  perfectly 
conducting  plates  with  separation  2H  located  in  a homogeneous  and  isotropic 
medium.  A Cartesian  coordinate  system  with  its  y-axis  normal  to  the  plates 
is  erected.  Both  plates  are  infinite  in  the  z-direction  and  finite  in  the 
x-direction  with  length  2L  as  shown  in  Figure  1.  All  figures  appear  at  the 
end  of  Chapter  II.  The  current  is  oriented  in  the  y-direction  and  is 
defined  as 

A 

J - y<5(x)  exp  (iSz)  , (!) 

where  3 is  the  propagation  constant  in  the  z-direction,  and  5(x)  is  the 
delca  function.  In  { 1 ] using  the  vector-potential  approach  and  the  Wiener- 
Hopf  technique,  we  obtained  a solution  for  the  problem  at  hand  in  a general 
form  for  any  parameters  with  one  restriction:  kL  must  be  much  greater 
chan  1,  i.e. , 

kL  » 1 , (2) 


where 

n T 

k * /a  eu  - 3 (3) 

and  £,  u are  the  homogeneous  media  parameters.  Using  the  solution  which 
was  obtained  in  (1],  we  will  perform  a numerical  calculation  for  the  case: 


W - | - 0. 16670 


(4) 


L 


3 


whe  re 


(5> 


Vb) 


o)/GU 


is  the  tree-space  wave  length.  3ecause  of  the  limitations  of  Equations 
(2),  (3),  and  (5)  we  calculated  numerical  results  for 

0 < 3 < 0.93  , (7) 

-a 

where  3 is  the  normalized  propagation  constant  5 - — — , It  can  be 

readilv  proved  from  Equations  (3),  1.4),  and  (3)  that  if  3 is  in  the  region 

( 0.93|,  where 

60  - 0.30008997  , (8) 

then  only  the  first  mode  can  propagate  in  the  x-direction;  therefore,  the 

electromagnetic  field  has  only  three  components  [E^,  H^_ , H_ ] . When  3 

intersects  the  point  3^  and  goes  to  zero,  the  second  mode  begins  to 

propagate,  and  the  EM-t'ield  consists  of  five  components  [E  , E , E . H , H 1. 

x y z x z 

Rewriting  the  tield  solution  from  [1]  for  the  case, when  no  more  than  two 
modes  can  propagate,  we  arrive  at  the  EM  field: 


E' 

(x,v ,z)  » 

E(x,y>  • exp 

( iSz) 

(.9) 

H’ 

( x , v , z ) - 

H(x,v)  • exp 

( i3z) 

(10) 

E(x,v)  - 

N N 

xE  + vE  ♦ , 
x y 

rE, 

(11) 
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xH  * zH 
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[M1+(k) )"  exp  (i2a) 


1 + -7=-  exp  (-  l -) 
v a -4 


(27) 


It  should  be  mentioned  that  we  investigated  the  lossless  medium  case; 
therefore,  in  the  region  < 3 < 0.93,  the  propagation  constant  for 
the  second  mode  has  onlv  an  imaginary  part.  Because  we  neglect  terms 
which  decrease  exponentially,  our  results  for  the  3bove  mentioned  region 
of  3 reduce  to: 


1)  E * E - 0 and 

x z 

2)  more  simple  expressions  for  the  other  three  components  of  the  field. 
We  3pply  numerical  analysis  only  over  the  regions  0<v<H,0<x<L. 


For  the  remainder  of 

correlations : 

the 

waveguide 

, one  can 

obtain  results  using 

the 

E (x,y) 

X 

* — 

Ex(-x 

,y) 

; E^Cx.y) 

* - E (x,-v) 

Eyu,y) 

3 

Ev(-x 

,y) 

; Ey(x,y) 

* E (x,-  ) 

E,(x,y) 

3 

E,  (— x 

,y) 

: E=(x,y) 

» - E^(x,-y) 

Hx(x,y) 

- 

H (— x 

X 

,y) 

! Hx(x,v) 

- H ( x , -v ) 

X 

H (x,y) 

2 

3 — 

H i-x 

,y) 

; Hz(x,y) 

= H,(x.-yl 

(28) 

It  is  interesting  to 

note 

that 

Ex: 

H are  continuous  and  that  E • 
2 y 

E , H 
z x 

are  discontinuous  when  3 crosses  3^  tor  more  exactly:  they  are  exponent iallv 
decreasing).  It  is  also  of  interest  to  determine  the  character  of  the 
behavior  of  the  x-component  of  Povnting's  vector.  As  one  can  easilv  see 
from  the  previous  expressions  for  the  EM  fields,  the  x-component  of 
Poyncing's  vector  for  the  second  mode  is  proportional  to  3^  and  goes  tv'' 
cero  when  3 -*■  3^.  From  this,  one  finds  that  the  energy  flow  in  the 

x-direction  is  continuous  when  3 intersects  the  point  3.. 

0 


6 
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Figure  5.  The  real  and  imaginary  parts  of  an  x-component  of  the  electric  field  as  functions 
of  a longitudinal  propagation  constant  for  points  of  view:  ^ = 0.9;  — = 0.0,  0 


figure  6.  ihe  real  and  imaginary  parts  of  a y-component  of  llie  electric  field  as 
of  a longitudinal  propagation  constant  for  points  of  view:  — = 0.1; 


Figure  H.  The  real  and  imaginary  parts  of  a y-component  of  the  electric  field  as^f unctions 
of  a longitudinal  propagation  constant  for  points  of  view:  — = 0.6;  — = 0.0,  0 
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;<u<l  imaginary  parts  of  a z-coinponent  of  the  electric  field  asyf unct ions 
gitudin.il  propagation  constant  for  points  of  view:  * = 0.4;  — = o.O  0 
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magi nary  parts  of  an  x-component  of  the  magnet  H:  field  as^f unctions 
n. 1 1 propagation  cons taut  for  points  of  view:  - = 0.6;  - = 0.0,  0 


The  real  anil  imaginary  parts  of  a z-component  of  the  magnetic  field  aSyfunctions 
of  .1  longitudinal  propagation  constant  for  points  of  view:  — = 0.4;  — = 0.0,  0 
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second  mode  is  responsible  for  this  misbehavior.  When  *r  * 0.3,  the 

H 

V 

contribution  of  that  step  decreases,  and  t'or  * 0.5,  it  equals  zero. 

One  can  observe  that  the  curves  are  smooth  it  3 = 3(,.  It  should  be 
mentioned  that: 

(. a)  The  discontinuous  behavior  is  observed  because  we  have 
neglected  the  exponentially  decreasing  terms; 

(b)  In  a lossy  medium,  the  rate  of  decrease  would  be  less.  The 
H -component  (see  Figures  18-11,  Ibl  doesn't  have  the  step  behavior 
because  the  contribution  of  the  second  mode  is  proportional  to  ct. , 
which  goes  to  0,  when  3 -*•  3^.  The  other  two  components , E^  and  E , are 
equal  to  zero  on  the  x-axis  and  have  their  largest  contribution  for 

y 

*-  = 0.5  (see  Figures  1-5,  10-13).  E,  is  smooth  and  E exhibits  Che  seep 

K X z 

behavior  at  3q.  The  figures  demonstrate  that  the  complexity  of  Che 
curves  occurs  approximately  in  the  region  0.7  < 3 < 0.8,  where  we 
observe  a sharp  peak,  which  is  due  to  the  resonance.  In  the  region 
0 < 3 < 0.5  the  figures  demonstrate  the  very  smooth  character  of  the 
curves . 

3 • Real  and  Imaginary  P arts  of  the  Field  Components  as  Functions 
of  Transverse  Coordinates 

In  this  section  we  present  the  real  and  imaginary  parts  of  the 
dominant  field  components  E^,  H , as  functions  of  the  x-coordinate 
for  two  values  of  ry  = 0.0;  0.5  and  two  values  of  3 = O.a;  0.~5.  The 
graphic  output,  shown  in  Figures  13-33,  was  obtained  using  the  results 
of  calculations  for  31  points  of  ^ in  the  region  [0.1,  0.A]  (step  = O.Oli 
All  Che  graphs  have  very  smooth  characteristics.  For  3 * A. '5  they  have 
slightly  more  complicated  form  than  for  •'  » ).«.  As  mentioned  in  the 
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previous  section,  for  “ = 0.5  the 

zero.  We  see  that  the  amplitudes 

v 

entire  region  of  view.  For  the  ^ 
the  eont ribut ions  of  two  modes, 
of  the  curves  are  changing  along 


contribution  of  the  second  mode  equals 
of  the  curves  are  constant  tor  the 
=0.0,  the  field  components  are  sums  of 
One  can  observe  that  the  amplitudes 
the  x-direction. 


I 


^1 


i 


xt< 


APPENDIX 


SOURCE  EXCITATION  OF  AN  OPEN,  PARALLEL-PLATE  WAVEGUIDE  PROGRAM 


'i 


t 


A complete  program  for  source  excitation  of  an  open,  parallel-plate  wave- 
guide program  is  presented.  The  computer  program  provides  three- 
dimensional  data-storage  for  the  real  and  imaginary  parts  of  five  components 
of  the  EM  field.  Data  were  obtained  for  ^ between  0.1  - 0.9  with  step  0.1; 

between  0.0  - 0.9  with  step  0.1;  and  3 - propagation  constant  between 
0.9  - 0.93  with  step  0.005  plus  (•)  0.30009.  These  data  were  used  to 
plot  EM  field  components  as  functions  of  the  propagation  constant.  The 
program  can  be  readily  modified  to  obtain  data  for  plotting  the  EM-field 
component  as  a function  of  the  x-coordinate . 
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r vnti.V  A.-i  At  I ELIU  INPUT. OUTPUT.  T APt:3  » T APE  1 --INPUT) 

DMF  ' LX  .» Bfv  « 3 UAL  * " l • F 1 » i I!  < % 0 » E Y 0 * E 2 0 * H X 0 . A 2 0 * C ON  S T » A A * S B 

: ■ : .'ll  NS  I ON  XX  V 9 ) . Y Y < 10 ) . 3ETAB  < 188'  - REXO  ( t0.9,l  J8  > . AMEXO ( i 0 
FKi  t > . 1 1)  . 9 » 1 OS  ' t AML  V .)  v 1 0 » 9 • 1 38  ' * KL  20  ^ 1 0 • 9 . L 88  ) « AMt  20  \ I 0 t 9 • 
i i\ u X 0 \ l 0 « 9 . l J !:! ' • AMHXC  v 10*9.  1 ! ! 3 ) » i < H 2 0 ( 1 0 » 7 « 138)  » AMHZO  C 10.9  * 
READ  t 1 . 2 ) W . T BK  . Y BPS  . r BRF  . XBR  1 , XBRS  » XBRF 
FORMAT  ^F7.5»«3<F7.3)‘) 

M = 1 000 

PI  = 3. 141592654 
CONST=CMPLXvO. . 1 . ) 

CON-2.  L 1593152 
BE  TAB ( l ) =0 . 

DO  X)  1 = 1.160 
BETAS  I If  l ) -BE TAB  v I > 10.005 
’0  CONTINUE 

BE  TAB v 162) =0.80009 
BETAS v 1*3 '=0.805 
BO  75  1=1.25 

75  BET ABC  I f 163)=BETAB( I f 162) +0.005 
CON T I NUE 
DO  98  N= 1.161 
BETA=BETAB(N) 

.\  = N 

A=10.*PI*5QRT< 1 .-B£TA**2) 
s=A*y. p [ 

AL.FAl=30RT  C B**2-  l . ) 

DD=PI*ALFAi  U 

CALL  BE 1 ^A.B.BBN.PI > M < CON • ALFA1 > 

CALI.  BE  2 < A * B . BBAL  » PI » M » CON  . ALF A 1 ) 

CALI  TE1 (A.B.T1 .BBK.PI ) 

CALL  .■  L2CA.B.T1.  ALFA  1 > DDK  * BBAL  * :■  1 < F2  ' 

X-XBR  I 
J=l 

30  Y=YBRI 

[=  t 

20  £XO=w;,[  *ALFA1  v 10.  *W**2'  ' *F2*SIN  v P I *Y  ' dSINvDI  *X> 

EYO=  10  . +P  I *<.  F l*COSC  A*X>  f < 1 . - C 1 . / ( 100  . *W**2  ' ) ) *F2*C0S  CP  I K Y 
iCOSaiB*X)  -CEXP^CMPLXCO.  . A*X>  ) ■ (2.  *A>  ' 

EZO=  -CONST  K BET  A*  \ PI  J ' * F 2 * 3 1 N C P I *Y  '.  *CQS  c DD*X  ' 

HXO=  - BETA*  10 . *F‘I  *<  FI  *COS  l A*‘<  > +F2*C0S<  Pl*Y  >*CQS<  j D * X ' - 1 l . 
iCEXPvCMPLXCO. . A * X ) ) > 

HZ0  = C0N3T  *A*F  1 *3  [ N v A*X  ) f CONS T KBB*F2  dCOS  vF  I K'  *8  I N «.  DD*\  ’ - 
fcO  . 5*CEXP  < CMF'LX  <.  0 . .A*X)  ) 

PE  XO  < I . J » i\  ' = R£AL  v EXO  ) 

AMEXO v I . J »K ) = A I MAG ( EXO ) 

RE YO  v I . J . is  ) =REAL  «,  EY  0 ' 


.9.138'. 
138 ) . 
133) 


THIS  PACE  IS  BEST  QUALITY  FRACTICABU 
JftQl  OOP!  IVBUUSHED  TO  DOC  -• 


ArfEYO  < I.  ,J»ls>  *AIMAG  < EYO ) 

REZO  s I » J.IS  ) -REAL  «,  EZO ) 

AMEZO  < I * J . Is  ) =AI MAG ( EZO ) 

RHXO  < 1 1 J » K ) -REAL ( HXO  > 

A M H X 0 v I • J . K ) =AIiMAG ( HXO ) 

RHZO< I » J » K ) =REAL ( HZO ) 

Ahl-iZO  i I * J • K ) -A  I MAG  ( HZO  5 
Y=Y  + YBRS 
1*1  + 1 

IF(Y.LE.YBRF)  CO  TO  20 

X=X+XBRS 

J=J  + 1 

TF(X.LE.XBRF)  GO  TO  30 
93  CONTINUE 

BO  99  L = 162 » 133 
8ETA=BETAB(L) 

K=L 

H=lO.*PI*SaRT<l . -BET  A#*2 ) 

X=XBR I 
J=1 

B=A*W/PI 

CALL  BE  < A » 3 . 3B  . PI .«  M ) 

X = XBR  I 

1 7 AA=CEXP  ( CMF'LX  (0,  . A#X  ) ) 

EYO*  -0 . 5*  < AA-BBACOS  s A*X  > )/SQRT( i . -8£TA**2) 

HXO=-BETA*EYO 

HZ0*O . 5*  < -SIGN  <1 . . X ) *AA+CONST*BB*S INC  A AX ) ) 

DO  LG  1*1.10 
REXO ( I » J » N ) =0 . 
mMEXO  ^ I » J . N ) *0 > 

REZO  C I . J . X ) =0 . 

AMEZO ( I » J * K ) =0 . 

REYC  < l » J » K ) *REAL  < EYO ) 

AME  (0  < I . J . is  ) * A I MAG  < EYO  ) 

RHXO  ( I * J » K > =REAL  HXO  ) 

AMHXOC I .J.  Is)  *A  I MAG;  HXO) 

RHZO ( I . J . N ) *REAL ( HZO ) 

1 3 a n H Z 0 < I » J « Is  ) =AIMAG  ( HZO ) 

X = X + X B R S 
J*J+1 

IFsX.LE.XBRF)  GO  TO  17 
99  CONTINUE 

UR  C TE(3»101)  REXO . AMEXO . REY  0 . AME YO . REZO  » AME  20  » RHXO . AMHXO . RHZO • AM 
iHZO 

101  FORMAT ( 10F3.3) 

0Y=0. 1 
DX*0. I 
TYs 1 ) =0 . 

DO  50  1=1.9 
YYs I + 1 ) =YY  < I ) +BY 
TO  CONTINUE 
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xx<.  t >=o.  l 

DO  *0  [=1.8 
xx v r t-i  > =xxv  i ) *-dx 
60  CONTINUE 

WRITE(o.lOL)  YY.XX.BETAB 
STOP 
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END 

3UBR0U TINE  BE 1 ( A . B . BBK . P I . M . CON  - ALFA 1 ) 

COMPLEX  F . BBK 
AML  =0. 

DO  10  I = 2 . M 

AM l =AM L + B/ I -AS I N v B/ I ) 

LO  CONTINUE 

F=CEXP(CMPLX<-B#PI/2. . B* < CON-ALOO ( B ) )-PI/2. fA+AMl ) ) 

BBN*< ALFA1+B>*F 

RETURN 

END 

SUBROUTINE  BE2 ( A » B . BBAL . PI . M . CON • ALFA l ) 

COMPLEX  F » BBAL • D 
AM  L =0 . 

DO  10  1=2  * M 

AM  l -=AM  L f ALF A l / I -AS  I N t.  ALFA  1 SORT  ( I **2-1  . ) ) 

LO  CONTINUE 

F=CEXP  ( CMPLX  ( - ALFA  1 *P  1/  2.  . ALFA  i *<  CON-ALOO  L B ) ) f )M i L-ALF A 1 *A/  S)  ) 
D=CMPLX ( L . . -ALFA1  ) 

BBAL  = F*D*S0RTL2. ) KALFA1,  B 

RETURN 

END 

SUBROUTINE  TE  t < A . B . T 1 . BBK .PI) 

COMPLEX  TL.BBK.Di 
D=B*SGR  H PI ) /SORT  v A*2 . > 

D L=CMPLX < 1 . f D » -D ) 

T t=BBK  K*2*Bl 

RETURN 

END 

SUBROUTINE  F12(A.B. Ti .ALFA l .BBK. BBAL. FI .F2> 

C 0 M P Li. X T l . BBK  » BBAL  * F .1  • F 2 . D 1 

01®  ALF  A L *<1  . fTl ) * ( BBAL  * *2  * t 4 . MIKS-  < t . «-Ti  )-<  ALFA  L f B ) X*2,  ALFA  L ' 
ii  2.  *ALFA  l )-L  . ) 

F 1 = i BBK  Mi2.  < C l . FT1 ) *A) ' ¥( 1 . f2 . #B¥&BAL*  K2  D L 
1-2=2.  *B*BBK.*BBAL/(A*Di  ) 

RETURN 

END 

SUBROUTINE  BE  t A . B » BB  * PI » M ) 

COMPLEX  T.F.BB 
CONST  = 1 . L 159313 16 
A M L = 0 . 
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In  Elect romagnetics  Laboratory  Report  No.  78-4,  we  nave 
calculated  the  electromagnetic  fields  as  functions  of  1,  the 
normalized  propagation  constant  in  the  z-direction.  The  para- 
meters chosen  for  the  computation  in  the  report  were: 


- » 0. 16670,  m = 5, 

L lo 


where 


w/cu 


At  the  request  of  Dr.  D.  Giri  of  SAI,  we  have  now  derived 
additional  numerical  results  for  the  following  choice  of  parameters, 
which  correspond  to  those  of  the  experimental  parallel-plate  structure 
being  investigated  at  Harvard. 


L = 12.5  tn 
H = 12.75  m 
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f = 25  Mny  (\^  = 12m) 


H 

L 


1.020  : 


1.04166' 
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The  propagation  constant  in  the  x-direction  can  be  written  in  the  form: 


a 
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1 isM. 

H T 


2 

- m 


1/2 

, m = 0,1,2, ..  . 


where  k = k^  / 1 - 3“ 


k = ^ 

0 :'o 


kL  <_  kQL  = 2tt  -r  = 6.344985. 

The  asymptotic  analysis  presented  in  our  report  78-4  was  based  on  the 
assumption  (kL>>!) . Consequently,  great  care  should  be  exercised  when 
the  range  of  application  of  these  formulas  is  extended  below  kL  = 10. 

It  is  not  difficult  to  prove  that  in  the  range  0 < 3 < 0.337916, 
only  three  modes  are  above  cut-off  in  the  x-direction.  Furthermore, 
two  modes  are  propagating  in  the  range  0.337916  • 3 0.382353  and  only 

one  mode  can  propagate  in  the  range  3 > 0.882353.  The  application  of 
the  formulas  and  computer  programs  developed  in  Report  No.  73-4,  though 
not  the  theory  itself,  is  restricted  to  the  range  where  two  modes  can 
propagate  in  the  x-direction.  For  this  reason,  we  develop  the  numerical 
results  only  for  the  region  3 _>  0.34,  and  specifically  for  the  range 

0.34  £ 3 <_  0.9. 

We  would  like  to  mention  that  it  is  possible  to  develop  the  necessary 
formulas  and  numerical  results  for  the  region  0 < 3 < 0.34  using  the 
theory  given  in  the  Electromagnetics  Report  No.  77-19. 
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t-igure  2.  Kea I and  imaginary  parts  of  E as  functions  of  B lor  parameters  given  in  the  lower  half  of 
Hage  1.  Ihe  y.,y  values  of  the  observation  point  arc  shown  in  the  inset. 
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I.  INTRODUCTION 

The  finite-width  parallel-plate  waveguide  is  used  as  an  EMP  simulator 
for  guiding  electromagnetic  waves.  There  have  been  a number  of  studies 
devoted  to  the  analysis  of  the  guide.  Rushdi  et  al.  [1]  and  Marin  [2],  [3] 
decompose  the  field  in  the  guide  into  two  components  (for  the  exp(ju,t) 
time  convention) : 

E(r)  = E1  + E2 

where 

\ = l An(x,y)e  J"n  , 

I2  - / /r  i(a,8)e-j(ax+Sy+Yz)dad6  . 

The  first  term  ^ describes  the  contribution  from  the  discrete  spectrum, 
i.e.,  the  so-called  "leaky-wave"  contribution.  The  second  term  is  a 
superposition  of  plane  waves,  the  contribution  from  the  continuous  spectrum. 
In  references  [1]  - [3],  only  E^  is  studied  in  detail.  This  information  is 
useful,  of  course,  only  if  E^  is  negligible  (for  certain  guide  geometries 
and  source  configurations).  To  test  this  assumption,  we  can  either  evaluate 
E2  directly  or  calculate  the  total  field  E instead.  We  take  the  second 
approach.  In  references  [4],  [5],  Krichevsky  and  Mittra  determine  E due 
to  a current  sheet  in  the  guide  by  the  Wiener-Hopf  technique.  Their  results, 
however,  are  valid  only  if  the  cross-section  of  the  guide  is  large  in  terms 
of  wavelength.  In  the  present  report,  we  again  calculate  E by  a different 
technique  moment  method  which  is  suitable  for  a guide  with  small-to-moderate 
cross-section. 


1 


for  different  d values,  and  obtain  the  total  field  by  superposition.  In 
this  report,  we  investigate  the  electromagnetic  fields  due  to  the  current 
source  in  the  form  of  J(x,y)e  J for  different  values  of  the  longitudinal 
propagation  constant  d.  The  step  involving  the  integration  over  the  spectral 
variable  6 is  not  carried  out  in  that  work. 
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II.  STATEMENT  OF  THE  PROBLEM 


The  geometry  of  the  problem  Is  shown  in  Fig.  1,  the  open  parallel 
waveguide  of  width  2L,  height  2H  is  infinitely  long  in  the  z-direction, 
and  is  excited  by  a y-direction  oriented  current  source  J1  defined  as: 


and  5 is  the  Dirac  delta  function.  Rect  is  the  rectangular  function  defined  as 

(l  i M £ 

rect  (t)  - < • (3) 

; elsewhere 

In  Equation  (1)  and  all  the  following  equations  in  this  report,  we  have 
adopted  a convention  in  which  all  equations  are  given  in  both  cases  of  odd 
and  even  mode  current  excitation  with  the  upper  equation  for  the  odd  case, 
the  lower,  the  even  case. 

The  total  electromagnetic  field  can  be  decomposed  into  two  parts: 
fC  - j?1  + ?s  , HC  » H1  + H8  . E1  and  H1  are  the  incident  fields  produced 
by  the  current  source  in  free  space  with  the  waveguide  removed.  Es  and  Hs 
are  the  scattered  fields  generated  by  the  induced  current  on  the  waveguide 
when  illuminated  by  E*  and  S1.  The  incident  field  is  first  examined  in  the 
next  section. 
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III.  INCIDENT  FIELD 

-M  * 

The  incident  fields  E and  H can  be  computed  via  the  v-oriented 

-<-i 

magnetic  vector  potential  A produced  by  the  impressed  current  source. 

fi  ' ,i/  x -J6z 
A - yA  (x,y)e  J 

The  z-direction  variation  for  all  fields  is  eJ  . Therefore,  it  is  suppressed 
in  the  subsequent  equations.  Also,  the  subscript  p in  a will  be  omitted  in 
the  following  equations  for  convenience  of  notation. 

The  vector  potential  A satisfies  the  Helmholtz  equation: 


32  . 32  . . 2]  ,i,  , Ti , . 

+ JjZ  + ^ A (x*y)  “ _uJ  (x*y) 


where  k2  * k2  - S2  ^ 0 , 


in  which  k » w/pe  is  the  wavenumber.  The  solution  to  (5)  is 
i-h  [sin  ay M 

A (x,y)  « 77  | ° < ) G(x,y  ; x ,y')dy’ 


-h  \cos  ay’ 


where  G is  the  Green's  function  defined  as 


G(x,y  ; x',y') 


(k®)  ; k-  > 32 


Jf  K0(kp  ; k=  < 8* 


and  s * /(x  - x')2  + (y  - y')2  . (9) 

/ \ 

and  Kg  are  zero-order  Hankel's  function  of  the  second  kind  and  modified 
Bessel's  function  of  the  second  kind,  respectively. 

and  {?■*■  can  be  derived  from  A*  by  using  the  following  equations 


ExU,y) 


. -i—  4-  *A  (x,y) 

juJUC  3x  Jy 


(10a) 
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Substituting  Equations  (7),  (11),  (12),  (14),  (15)  together  with 
Equations  tS)  and  1.13)  into  Equation  tlO),  we  can  obtain  E*  and  it*. 
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IV.  SCATTERED  FIELD 


the  Induced  current  J3  on  the  waveguide  has  two  components. 


■ts  * ,s  , . ■ ,s  . 

J ■ xj  u)  * iJ  u)  ; 

X 2 


^ Lt> ' 


the  magnetic  vector  potential  A&  due  to  the  induced  current  J3  also  has 


two  components. 


xA*(x,y)  + 2AS(x,v) 


(17) 


s s 

Both  A^  and  A,  satisfy  the  Helmholtz  equation: 


. s » \ 

, ~ “ + 77^  + AS(x,y)  - -uJSix)  , 
\,3x-  lv-  ti  u u 


(,  13) 


where  u ■ x or  2,  and  the  solutions  are 


rL 

u 

| 

J-L 

•L 

*J 

J-L 

and  H can  be  derived  from  \s  by  using  Che  equations  given  below 


s . . 1 3*  s 3 « 


E‘  vx.y>  - - J.oA'  tx. v>  - .. 

x x jotue  i 3x-  x 


V**y>  - -,s  Vx,y>i 


CU) 


E%Stx.v) 


lT  I TTw  A!vx.y)  - JS  ^ A^x.y't 


w3x3y  \ 


E3(x,y)  " -J-uvfix.v)  j ASix,y)  - JSA*(x,v)  i 

* a-  JolUt  ■ * X X a- 


. . s . 1 ) ,s  . 
HxiX.y)  - - — Aau.y> 


H^ix.y)  - ~ [-JdA3(x,y>  - A*tx,y)] 


Ulb) 

t:ic> 

vCld) 

(.  Sle) 


H ix.vt 


L 3 .s  . 
- — A ^x.v3 
u 3v  x 


I ] 
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electric  field  be  zero  on  the  surface  of  the  waveguide,  viz., 

E*(x,h)  + ES(x,h)  - 0 ; |x|  <_  L (22) 

E^(x,h)  + E^fx.h)  *0  ; |x|  <_  L (23) 

Due  to  the  symmetry  of  the  problem,  we  do  not  have  to  enforce  the  boundary 


condition  on  the  lower  plate  (y  * -h) . Equations  (22)  and  (23)  can  be 
explicitly  written  as 

+ k.^1  A^(x,h)  - jt3  — A*(x,h)  - -j,aueE*(x,h)  ; j x j <_  L , (24) 

(-3‘  + k-)A3(x,h)  - Jfi  ~ AS(x,h)  * -juiucEi(x,h)  ; |x|  < L (25) 

Z J X X 2 

We  then  substitute  Equations  (19)  and  (20)  into  Equations  (24)  and  (25)  and 

numerically  solve  for  J*’  and  J*  to  obtain  E and  H . 

X z s s 


It  should  be  pointed  out  that  the  A's  In  Equations  (28)  and  (32)  are  not 

necessarily  the  same  3S  those  in  Equations  (26)  and  (27);  we  have  chosen 

the  same  symbol  for  the  sake  of  convenience.  Next,  we  enforce  (26)  and  (27) 

simultaneously  at  a set  of  N match  points  (x  'sh  Again,  (x  's;  are  chosen 

m m 

to  coincide  with  those  points  ix  's;  in  (29)  for  the  sake  of  simolicitv. 

n 

We  then  have 


Ax(xm+rh)  * ‘ :)Ax(Vh)  + Ax(Vrh) 


-)  f |A*ivrh>  - A;<vi-h)|  ■ -i  ~ 4(\m 


(34) 


and  -j  [A^tx  ,h)  - AS(x  . ,h)]  + k:d:As(x  ,h) 
2 x tn+1  x m-L  t z m 


-j  — E*(x  h) 
jj  z m 


(35) 


where  m * 1,  2,  • • • , N. 

3y  substituting  Equations  (28)  and  (32),  the  expansion  function  expressions 


for  J®  and  j®,  into  Equations  (19)  and  (20),  AS(x  ,h)  and  AS(x  ,h)  become 
x - x m z m 


N 


As(x  , h)  - ~ Z,  a Sa 
x m 4j  n*l  n m, n 


(36) 


where  Sa 


m.n  I 


x +A/2 

[G(x  ,h  ; x',h)  ± G(x  ,h  ; x',-h)]dx'  ; 

x -A/2  m m 

n 


(37) 


n * 1,  2,  •••,  N, 
m = 0,  1,  • • • , vN+l) , 


s u N b 

and  A*(x  ,h)  - ~~  Z,  b Sb 

s m ■*.)  n*l  n m,n 


(38) 


b !x  i+A/2  “a — 

Where  Sm,l  * ‘ /'rT^T  [G(xm’h  ; x’’h)  1 G(xm’h  ; > (39a) 

* ™ la 


Sm.n  * S»,n  5 n » 2,  3,  • • • . (N-n  , 


( 39b) 
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(39c) 


SD  M = i — ~rf  [G(x  ,h  ; x',h)  ± G(x  ,h  ; x',-h)]dx' 

Jv  _A/2  ^ L “ X m 


V 


m = 0,  l,-*-,  (N+l)  . 

Substituting  Equations  (36)  through  (39)  into  Equations  (34)  and  (35),  we 
arrive  at 


[S*  + (k2A2  - 2 )Sa  + Sa  ][-j  ^ (Sb  - Sb  )] 

nrfl,n  m,n  m-l,n  2 m+l,n  m-l,n 


l-J  f <Cl.»  - 


[k2A2S°  ] 

t m,  n 


<bnJ 


^ Ai 


& Bi<vh)i 


[~j SB  (xm*h) ] 


(40) 


The  above  equation  can  be  readily  solved  for  a and  b^.  The  unknowns  A^  and 


A3  can  then  be  calculated  using  Equations  (36)  through  (39),  and  the  scattered 


field  can  be  determined  from  Equation  (21)  and,  therefore,  the  total  field. 


The  equations  for  calculating  ^'s  and  ^'s  are  given  in  the  following: 


(A)  Evaluation  of  S 's 
m,n 


(i)  If  m = n 


S = I„  ± image 
m,n  0 


(41) 


where  image 


fX  +A/2 

G(x  ,h  ; x' ,-h)dx' 

x -A/2  m 

n 


(42) 


rX  +A/2 

and  I-  = I m G(x  ,h  ; x',h)dx' 
0 Jx-A/2  m 


(43) 


m 


= 2 


x +A/2 


G(x  ,h  ; x’,h)dx' 


(44) 


- x )]dx’  5 > 


0 1 t 


;X  -hi/2  , , 

iti  j m K0[kt(x'  - xm)]dx'  ; k2  < ft- 


(45) 


-=-  [ C H^2)(t)dt  ; k2  > S* 

kt  ^0 


(46) 


rk  A/2 

— c K (t)dr.  ; k2  < B* 

nkt  Jo  ° 


Since  simple  algorithms  have  been  developed  Co  evaluate  the  integrals 
of  HU)  and  Kq  from  zero  to  any  positive  number  [6],  Equation  (46)  can  be 
evaluated  efficiently  and  accurately* 

(ii)  If  ® ^ n , 


r 


s3  - { 

m,n  \ 


x -A/2 
n 


G(x  ,h  ; x'  ,h)dx'  - I “ G(x  ,h  ; x\h)dx*  i image  ; m>n 
m Jx  +A/2 

n (47) 


rx  +A/2 

n G(x  , h ; x',h)dx’ 
I m 

* V 


;x  -A/ 2 

n G(x_,h  ; x’  ,h)dx’  i image  ; m<n 
x 

m 


m 


where  the  image  is  given  in  Equation  (■+-) 
(B)  Evaluation  of  Sb  's 

>a»  U 

First  consider  Sb  which  is  expressed  as: 
ni,  i 


(43) 


S i “ I,  + 12 
m,  1 1 


-L+A'2 


where  I 


1 J _L  v L + x m 


r (G(x  ,h  ; x’.h)  ± G(x  h ; x\-h)]dx’  <«) 


(X,+A/ 2 


1 


tj  - J^-  r , (c(x„,i.  . «•*»  = «vh  '•  x'--h,ldx'  - '50' 
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[*■  ' 


rw 


(i)  If  in  + 0 , 


[G(Xm,h  ; x',h)  t GCx^.h  > x ' , -h ) ] 


f-L«/2  r /r— 

L1  J_L  V'  L + x' 

- /rhr  [G(vh  : -L*h)  1 G<vh  : -L.-ioj 


dx' 


+ [G(x  ,h  ; -L,h)  ± G(x  ,h  ; -L,-h)] 
m m 


-L+A/2 


-L 


where 


( 


-L+A/2  r 


! f 1 ' djt' 

> L + x 

( 5i ) 


j / L + x1  dx'  Can  evaluaced  analytically, 


which  is 


■L+A/2 


/ A 


-L 


v L + x 


- dx'  * /2  A . 


(51a) 


(ii)  If  ra  - 0 


first,  the  small  argument  behavior  of  function  G is  obtained: 


_ ii 


lim  G(-L,h  ; x',h)  “ - — *+  Ln 
x'*-L  " 


Ykc(L  + x') 


(52) 


where  y ” 0.5772156649.  Then,  in  evaluating  1^  in  Equation  (49), 
the  singularity  at  x'  - -L  is  first  subtracted  from,  then  added 
to  the  integrand,  and  by  using  Equation  (52),  we  have 
-L+A/2 


-L+A/2  , ( 

,L  /rrr 

•>,  frt  (l  * x'llA 

f ‘■L-Hr- J j 


h)  i G ( - L , h ; x',-h) 


■dx' 


^ |_L  /rTTr("':n(L  + x’> + d> 


dx  ' 


( 5J) 


The  first  integral  in  the  above  equation  is  amenable  to  numerical 
integration,  and  the  second  integral  can  be  evaluated  analytically: 


15 


: i 


-L+A/2 


- ii 

IT 


/ A 

/ : — : — y tn 


/ L + x' 


rkc(L  + x') 


dx' 


¥ a 4 {”  - 2}  • 


- - 11  /i 

IT 

On  the  other  hand,  I,  can  be  evaluated  as 
I 


fX^A/2  r r—^ 

2 ■ i /mp-  ^G(Xm,h  ; x'*h)  ± G<X  -h  • x'.- 

Jx1-A/2  V L + x m m 

- [G(x  ,h  ; x'.h)  t G(x  ,h  ; x',-h))\dx'  + Sa  , 
ra  ra  J m,  1 


m ■ 0,  1,  • • • 1 (N+l) 


Substituting  the  above  expressions  for  1^  and  I,  into 


Equation  (48),  we  then  obtain  Sb  . Following  similar 

ni)  i 


procedures,  we  obtain  the  equations  for  Sb  : 

m,  N 


Sm,N  X3  + X4 


where  I, 


L (/^ 

L-A/2  \ L x 


C G C x , h ; x',h)  t G(x  ,h  ; x',- 


»'  L - x 


r (G(x  ,h  ; L,h)  ± G(x  ,h 


; L,-h)]J 


dx' 


+ [G(x  ,h  ; L,h)  ± G(x  ,h  ; L,-h) ]/2  ^ , 


in  ' mr 

when  m * 0,  1,  • • • , X 


or  I 


f 


3 J 


h) 


/ r~rr  (c(L,h  ; x’,h)  ± G(L,h  ; x’,- 
L-A/2  L x l, 

2j  , >kt(L  ~ x'^  . 21  - f PM*!  >| 

it  cn — 2 — Jdx  - '2  - 2j 


when  m ■ N + 1, 


In 


m - 0,  1,  •••  (N  + 1)  . 


Using  the  equations  obtained  above  to  compute  Sa  's  and  Sb  's  in 

m,n 

Equation  (40)  and  inversing  the  matrix,  we  then  are  able  to  determine  the 

a ' s and  b ' s . 
n n 
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VI.  NUMERICAL  RESULTS  AND  DISCUSSION 


1 

I 

1 

I 


The  results  evaluated  by  the  methods  discussed  in  the  last  section 

are  now  presented  in  graphical  form.  The  physical  dimensions  of  the 

waveguide  and  the  current  source  for  which  we  perform  the  computation 

are  (see  Fig.  1 and  Equations  (1)  and  (2)) 

L - 12.5  m 

H ■ 12.75  m 

h • 5 m 
o 

x • 0 m 
o 

? - 1 • 

The  frequency  used  is  25  MHz,  and  the  current  source  is  of  even  mode. 

The  magnitudes  of  the  induced-current  components,  j and  |3*| , are 
plotted  in  Figs.  3(a)  and  3(b)  as  functions  of  the  transverse  coordinate  of 
the  waveguide.  They  are  presented  for  different  3/k  values  in  the  range 


between  0.86  and  0.99. 

The  plots  of  the  magnitudes  of  the  electrical  field  components,  |£  I , 

J?  [,  and  |£  ; , consist  of  two  parts. 

(A)  Field  components  as  functions  of  3/k: 

E , E , and  E are  plotted  for  3/k  in  the  range  between  0 and 
x y 2 

0.999  in  Figs.  4(a),  4(b)  and  4(c),  when  x - 0.5L,  y » 0.47H. 

It  is  observed  from  these  plots  that  for  most  values  of  3/ k between 
0 and  1,  |Ex| , jE^J,  and  |E_J  are  relatively  small  and  vary  rather  slowly 
with  respect  to  3/k;  however,  for  3/k  in  the  region  between  0.7  and  1, 
the  magnitude  of  the  field  component  is  found  to  have  a sharp  peak. 
Therefore,  the  major  contributions  to  the  total  fields  should  come  from 
this  portion  of  the  spectrum. 


■ 
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0.8  -0.6 

-0.4 

o 

CM 

o 

0.2  0.4  0.6  08 

x/L 

Figure  3(a) 

. ij;i 

as  a function 

of  x coordinate. 

1.2  -1.0 


-0.8  -0.6  -0.4  -0.2  0 0.2  0.4  0.6  0.8  1.0 

x/L 


Figure  5(c). 


El  as  a function  of  x coordinate  for  y/H  = 0.47. 
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APPENDIX 


A computer  program  for  determining  the  magnitude  of  all  six  component 
of  the  EM  field  is  presented  here.  Data  were  obtained  for  x/L  between 
0-1  with  step  0.05,  3 between  0.379  - 0.999  with  step  0.02,  and  y at  any 
given  value.  The  program  can  be  readily  modified  to  obtain  the  real  and 
imaginary  parts  of  all  field  components  as  functions  of  the  x-coordinate 


and  3 . 
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